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Abstract 

Charmonium productions in p—p and A — B collisions have been estimated 
within the ambit of colour evaporation model (CEM). The model parameters 
have been fixed by fitting the theoretical results with CDF data. The method 
is then applied to RHIC and LHC energies to obtain the transverse momentum 
distributions of J/i(i, ip' and Xc- Suppression due to Debye screening in a 
quark gluon plasma (QGP) is estimated at various centrality cuts. The final 
Pt distributions of various resonances are then predicted convoluting with the 
survival probability. 

I. Introduction 

Ever since the possibility of creating quark gluon plasma (QGP) in relativistic heavy 
ion collision was envisaged, numerous signals were proposed to probe the properties 
of such an exotic state of matter. In this context Satz and Matsui [1] had suggested 
that the production of heavy quark resonances {J/ip) will be suppressed as a result of 
colour Debye screening in a hot and dense system of quarks, anti-quarks and gluons. 
This suppression could be detected experimentally through the dileptonic decay 
mode of these resonances. ALICE dimuon spectrometer [2] is dedicated to look for 
this type of signal. However, it is a daunting task to disentangle the contributions of 
the heavy quarkonium states to muon spectrum due to the background from several 
other sources, e.g. Drell-Yan, semileptonic decay of open heavy flavoured mesons 
(DD, BB) etc. Low energy muons from kaons and pions also constitute a large 
background. 

In this work we shall estimate the hard charmonium productions {i.e. J/ip 
produced from initial hard process, will be called hard J/ if) hereafter) both in p — p 
and A — A collisions. The production of heavy resonances proceeds via the following 
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two steps: (i) the production of heavy quark- ant iquark pairs (perturbative), (ii) their 
resonance interactions to form the bound state (non-perturbative). 

The initial state in relativistic heavy ion collisions consists of either hadronic 
matter or QGP depending on the incident energies of the colliding nuclei. At LHC 
energies the formation of QGP is unavoidable. Even if the system is formed in QGP 
phase it will revert back to hadronic phase due to the cooling of the expanding 
QGP system and hence the interaction of the hard J/ijj with the hadronic matter 
is inevitable. Therefore, in addition to the suppression due to Debye screening, 
one needs to consider the survival probability of those J/ip due to its interactions 
with the hadronic medium. In this work we have neglected the suppressions due to 
co-movers as this effect is found to be negligibly small. 

The paper is organized as follows. In section II we shall describe the formalism 
for charmonium productions in the CEM along with the survival probability as a 
function of centrality. In section III results of our calculations will be presented 
followed by summary and discussion in section IV. 

II. Charmonium Productions 
a. Nucleon-Nucleon Collisions 

The CEM (also known as local duality approach) [3, 4, 5, 6, 7, 8] states that a heavy 
quark pair with mass Mqq < IMqq transforms (q designates light quark, while Q 
stands for c or b quarks ), independent to its colour and spin, to a QQ bound state. 
The production of a particular bound state depends on the dynamical details of 
hadronization procedure. The bound state formation probability (F[nJ PC ]\ the so 
called normalization factor) cannot be calculated from first principle. This is treated 
as a parameter and can be extracted by fitting the model with the experimental data. 
F[nJ PC ] depends on the particular resonance state under consideration, mass of the 
heavy flavour, the order (LO or NLO) of the QQ production and the type of parton 
distribution function (PDF) used. Another important point is that the produced 
QQ is not constrained to have the proper spin-parity and colour neutrality. The 
pair sheds its colour non-perturbatively to evolve into asymptotic state without 
affecting the cross-section. This model cannot predict the absolute cross-section of 
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heavy resonance production, it however, predicts, their and y/s dependence. As 
mentioned, CEM does not bother about the spin-parity and colour neutrality of the 
bound state. This constraint can be removed by considering an improved version 
of the duality approach, where one assumes that only the colour singlet part of the 
cross-section contributes to the bound state production. It might be recalled here 
that the CEM contains some features of non-relativistic QCD, due to the inclusion 
of colour octet processes. 

In spite of these limitations of CEM, it is capable of explaining p? distribution of 
charmonium in hadron-hadron collisions reasonably well. As mentioned before, the 
production of charmonium consists of two stages; production of a cc pair (pertur- 
bative process) and subsequent non-perturbative evolution into asymptotic states. 
We have considered those hard processes which can contribute to cc productions 
irrespective of their colour and spin-parity. The colour neutralization occurs by the 
interactions (one or more soft gluon emission) with the surrounding colour fields and 
this step is considered to be non-perturbative. In CEM quarkonium production is 
treated identically to open heavy flavour production with the exception that in the 
case of quarkonium, the invariant mass of the heavy quark pair is restricted below 
the open charm/bottom mesons threshold (see eq. (2) below), which is twice the 
mass of the lowest meson mass that can be formed with the heavy quark. Depend- 
ing on the quantum numbers of the QQ pair different matrix elements are needed 
for various resonances. The effects of these non-perturbative matrix elements are 
combined into the factor F[nJ PC ] which is universal [9] i. e. process and kinematics 
independent. It describes the probability that the QQ pair forms a quarkonium of 
given spin (J), parity (P) and charge conjugation (C). The production cross section 
for a J/j/i or any charmonium state is therefore given by [9] 



where the non-perturbative (long distance) factor can be written in terms of the 
probability to have colour singlet state (1/9) and the fraction p R of each specific 
charmonium state. The perturbative contribution (short distance) is given by 



a[R{nJ PC )} = F[nJ PC ] d[QQ] 



(1) 
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The contributions to heavy quark production in leading order come from qq — > 
QQ and g g — > QQ. The differential cross-section for heavy quark pair production 
in hadron-hadron collision is given by [fO] 

[h A h B -> QQX\ = , (3) 



dM 2 Q dy 



where 



a qq^QQ 



H(x a ,x b ,Q 2 ) = q h f A (xa,Q 2 )q h f B (x b .Q 2 )+ q h f A (x a , Q 2 )q h f B (x b , Q 2 ) 

f 1 J 

+g hA (x a , Q 2 )g hB (x b , Q 2 )v gg ^ Q Q- (4) 



x a>b = Mqq e ±y j y/s, y/s being the centre of mass energy of the hadronic system, 
y stands for rapidity and Mqq is the invariant mass of the pair, q/s and g's are 
the PDFs for quarks and gluons respectively, these are to be taken from either 
CTEQ or MRST or GRV [11]. Results presented in this work have been obtained 
with CTEQ(LO) distribution function. Combining eqs.(l),(2) and (3) we obtain the 
cross-section for resonance production per unit rapidity as, 

da n r, r 2m D/B „ H(x a , x b , M^n) 

-[h A h B ^ RX] = F[nJ PC ] j ' dM 2 QQ 1 ^ . (5) 

The above equation can be used to calculate the longitudinal momentum de- 
pendence of the quarkonium production cross section as shown in Ref. [10]. 

The leading order (LO) calculation does not have any dependence as the 
heavy quark pairs are produced with px = (assuming there is no px broadening 
of the initial state partons). In order to obtain the px distribution one has to go 
beyond LO. The dominant production mechanism of heavy quark pairs with large 
Pt and invariant mass near the threshold is the large pt gluon splitting with the 
probability given by [5] 

dProb a s 1 

dM^=6^Mf Q () 

and 



d 2 p T dy hm,Q QQ d 2 p T dy 6tt M 2 n ' 
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where N stands for nucleoli, da 9 / cPprdy is the inclusive gluon px and y distribution 
in NN collisions calculated in LO. It is given by (a b — > gc) 

J^r y -Te^E / v( I .,« ! )G 1 ,»( Ii ,e ! )«.g, (8) 

where x a = pr(e y ' + e y4 ) / 'y/s, Xb = Pt{.^~ % ' +e~ yi ) j ' \fs, (Ai) 2 is the matrix element for 
the process and Gs are the parton distribution function, px and y are the transverse 
momentum and rapidity of the gluon which splits into QQ pair and 2/4 is the rapidity 
of the particle c. The processes that contribute to the heavy quark pair productions 
are qq — > g g, g g — > and g q(q) — > gq(q)- Here the final state gluon splits into 
a heavy pair. Our analyses show that F[J/ip] = 0.045 reproduces the CDF [12] data 
quite well (see later). 

b. Nucleus-Nucleus collisions 

The charmonium production, is observed to be suppressed both in p-A and A-B 
collisions (compared to the scaled p-p scattering). This could be either due to nuclear 
absorption or due to the reduction of Q-Q interaction range in a QCD plasma caused 
by Debye screening. While the former is known as the normal nuclear suppression, 
the latter driven by the plasma effect is dubbed as anomalous suppression. 

To understand the mechanism of anomalous suppression one introduces the 
concept of quarkonium formation time (r /) and the dissociation temperature Td 
determined from the condition at which the Q-Q interaction range becomes equal 
to the size of the quarkonium. The corresponding time when plasma attains a 
temperature T = Td, is denoted as Td- r / on the other hand is the time required 
for the Q-Q pair to evolve into a physical charmonium state. This in the plasma 
rest frame would be Lorentz dilated and reads as r form = 7T / = r f\Jl + p T /m 2 R . 
High pt quarkonium states can evade suppression under two circumstances: (i) 
if Q — Q pair materializes into a bound state when the plasma has cooled down 
below Td or (ii) when the bound state is formed outside the plasma. The first 
condition, i. e. r form > r d implies no suppression for p T > pf^ while px > P2T 
for no suppression stems from the condition, \ f+ T fp T /rrifi\ > R s . Here R s is the 
radius of the screening zone (see later) and r is the position where the heavy quark 
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pair is produced. Therefore anomalous suppresion will be realized for quarkonium 
momenta less than mm{p^p, P2t}- 

Next we consider J/ijjfy', Xc) production in p — A and A — B collisions. To this 
end, we first briefly mention the necessary formulae in Glauber model [13, 14]. The 
total inelastic cross-section in A — B collisions at an impact parameter b is given by 



da 



AB 
in 



1 - 



1 - ^(b) <7j 



NN 



AB 



= 1-P (b), (9) 



dh 

where T A B(h) is the nuclear overlap function given by 

T AB {h) = J dsT A (h)T B (h-s) (10) 

The nuclear thickness functions are normalized to unity, i. e. J dhT A (h) = 
J dh dz p A (h, z) = 1, here p A is the nuclear density distribution. 

Generally we are interested in the cross-sections for a set of events in a given 
centrality range defined by the trigger settings. Centrality selection corresponds to 
a cut on the impact parameter, b of the collisions. The sample of events in a given 
centrality range < b < b m , contains a fraction of the total inelastic cross-section. 
This fraction is defined by [15] 

Now we discuss the J/ip survival probability when, after production, it prop- 
agates through the target /projectile nucleus. After creation, the J ftp meson can 
interact with other nucleons in the target and the projectile and may get destroyed 
mainly due to J ftp — N interactions. The cross-section for J/ip production in p — A 
collisions can be written as [16] 



a j/^( b m) = ^ J o dhdzp A (h,z) exp —(A — I) J <j abs p A (b, z') dz' 



(12) 



where a™ is obtained from eq.(5). The interpretation of the above equation is 
as follows. The resonance is formed at f = (b, z) where the density of the target 
nucleus is Pa(j)- It can travel in forward direction (z) at constant impact parameter 
and its intensity is attenuated due to J/ip — N inelastic collisions. The exponential 
factor accounts for this attenuation. 
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The generalization of eq.(12) in nucleus-nucleus collisions is straightforward. 
The J / ip production cross-section in A — B collisions at an impact parameter b can 
be written as 



The observed charmonium suppression data at SPS energies [17] could be ex- 
plained either by the Debye screening in a QGP or by co-mover scattering due 
to hadronic matter in the initial state [16, 18]. But the co-mover interpretation 
of heavy quarkonium suppression at LHC energies does not seem to be plausible, 
as the initial temperatures achieved at these energies could be considerably larger 
than the transition temperature, T c . Thus, the heavy quarkonium will encounter the 
hadronic comover at a much later time. It should be noted here that the formation 
of QGP does not necessarily guarantee the suppression of the heavy resonances. 
The temperature of the system must be greater than the dissociation temperature 
of the particular resonance in order to get suppressed. The dissociation temperature 
has recently been calculated using lattice QCD [19, 20]. The most recent values of 
T d for different heavy flavour resonances are shown in Table I. Thus as long as the 
plasma temperature remains greater than T d bound state cannot be formed. The 
time Td is obtained by Bjorken scaling law [21]: 



where Tj and Tj(6) (defined later) are the plasma formation time and initial tem- 
perature of the plasma respectively. In order to obtain the centrality dependence 
of the suppression we estimate the initial temperature by assuming the isentropic 
expansion of the system, namely, 




(13) 



Td(6) = Ti 



Ti(b) 



(14) 
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Resonance 


m 1 m r-i 

T d /T c [19] 


m 1 m r-i 

T d /T c [19] 


(T d )/T c 


T / / 

J/ip 


1.1 


2.0 


1.5 


Xc 


0.74 


1.1 


0.9 


ip' 


0.1-0.2 


1.1 


0.625 


T 


2.31 


4.5 


3.4 


Xb 


1.13 


2.0 


1.55 




1.1 


2.0 


1.55 


Xb 


0.83 








0.74 







Table I: Dissociation temperatures of charmonium and bottomonium system in pure 
gluonic plasma. 



T?{b) 



2n 4 



dN AB 



(b). 



(15) 



45C(3) Aa^R^n dy 
a,i = 7r 2 /90 (21AT F /2 + 16), is the degeneracy in the initial state, dNAB/dy is the total 
multiplicity of produced hadrons measured experimentally and R t is the transverse 
dimension of the plasma.. One would expect that the total multiplicity in nucleus- 
nucleus collisions at RHIC and LHC energies will come from both soft as well as hard 
collisions. It has been shown recently [22] that about 10 % of the total multiplicity at 
RHIC energies comes from hard collisions. Therefore, the total hadron multiplicity 
in nucleus-nucleus collisions at a given impact parameter can be written as (the so 
called two-component model) 



dNAB 

dy 



(b) 



[1-x) 



iV part (6) 



+ xN coll (b) 



dN, 



NN 



dy 



(16) 



In the above x is the fraction of hard collisions. For hadron rapidity density 
(dN NN /dy) in nucleon-nucleon collision we use the following form [23]: dN NN /dy = 
1.5(2.5-0.25 ln(s) +0.023 ln 2 (s)). iVp art (6) and N coll (b) are the numbers of partici- 
pants and number of collisions respectively at an impact parameter b. In non-central 
collisions the plasma transverse radius can be approximated as Rt ~ 1.2[0.5A^ part (6)] 1 / 3 [24]. 
The average initial temperatures and the dissociation times for various centrality 
classes have been shown in Table II. 

The critical radius R s , beyond which there will be no suppression, is given by [25] 

1/2 



Re — Rt 



(17) 
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Centrality (fm) 


Rt (fm) 


(Ti) (GeV) 


fa) (M 




[0-3] 


6.9 


0.617 


4.43 




[3-6] 


6.36 


0.582 


3.73 




[6-9] 


5.38 


0.522 


2.70 




[9 - 12] 


4.0 


0.438 


1.59 





Table II: Initial temperatures, transverse radii, and dissociation times at different 
centralities for J ftp with T d ~ 1.1T C and number of flavour, N F = 0. All the averages 
here are over the impact parameter. 

The survival probability, which is the ratio of bound states produced by the 
escaped pairs relative to the number of bound states that would be formed in the 
absence of QGP, is given by [25] 

3 C 0m 



> - (*)■ 



1/2 



r dr 



(18) 



where 



0r 



= 7r : 
= cos~ 
= : 



z < -1, 
- 1 < z < 1, 
z > 1 



(19) 



with z = [(R 2 — r 2 )mn — r^p\ /ttir]/ {2rro/Pr) • We assume a pure gluonic plasma 
with T c ~ 270 MeV [26]. 

Suppression due to Debye screening continues till the temperature of the QGP 
drops below T d . The QGP starts hadronizing at T c (equivalently, at tq = Tfri/T^, 
the time when phase transition starts). During the time interval tq — r<j, the re- 
maining resonances will not be dissociated as long as T c is lower than T d . During 
the mixed phase interval (th — tq) QGP part does not contribute to the suppression 
for T c < T d . Here th = QqTqI QH-, where qq and qh are the statistical degeneracy for 
QGP and hadronic phases respectively. However, the co-moving absorption starts 
at tq and the density of the co-moving hadrons is determined by the temperature 
T c . Following Ref. [27] we have checked that this effect is very small in the present 
scenario. Hence we do not consider suppression due to co-mover absorption here. 

Finally, the pr-distribution of heavy resonances is given by 
dNj/^ _ da^f^ d%S Q (p T ,b)ABT AB (b) 



dp T dp T / fem d 2 b [1 - P (b)] 
9 



(20) 
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Figure 1: CDF data [12] for prompt J/i/j production is compared with CEM (LO 
plus gluon splitting) prediction at yfs = 1800 GeV. CTEQ(LO) has been used in 
all the calculations. 



However, with nuclear absorption the above equation will be modified to: 

dNj/ip 



rlrr AB 

Io m ^d%S Q ( PT ,b) 



dp T j° m d% [1 - P (b)] 

Similar expression can be obtained for if)' and Xc- 



(21) 



III. Results 

In Fig. (1) we compare our results with CDF data where m c = 1.3 GeV, Q 2 = m^+p 2 -, 
have been used. It is clear that the data is reproduced reasonably well with the values 
of F given in Table III. Any effect due to /cy-broadening is ignored here. However, 
it might be mentioned that the low pt data on T production could be explained in 
color evaporation model by assuming substantial kr smearing [7]. 



F[J/if,] 


F 




Ej B( X cj - J/i)X)F[ Xc j\ 




0.045 


0.0126 


0.005 





Table III: Values of the universal factor used in the calculation (no feed down). 

Next we consider the J /if) production at LHC energies. First of all, in Fig. (2), 
we show da/dp? of J/if> and if)' from p — p collisions at ^fs = 5500 GeV. Here the 
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io'ft 




P T (GeV) 

Figure 2: Differential cross-sections of J/ip and ip' in p — p collisions at y/s = 5500 
GeV. Only direct contributions are shown. 

yields correspond to the direct productions i.e. feed down from higher resonance 
states have been ignored with m c = 1.3 GeV. The feed down contribution can be 
estimated by multiplying the yield with appropriate branching ratio. 

The differential pt distributions of J/ ip and ip' from p — p scattering (as shown 
in Fig. (2)) provide the baseline for Pb — Pb collisions at LHC energies. To obtain 
charmonium pr distributions at LHC energies we use eq. (20) which gives pt spectra 
without nuclear and comover absorptions. The amount of nuclear absorption can be 
calculated using a abs = c r abs ( v / so) {s/s ) A/2 , where s = 17.3 GeV, 0- abs (s o ) = 5 ±0.5 
mb and A = 0.125 [28]. The nuclear absorption of heavy resonances at SPS energies 
is well studied and at LHC energies it can be eliminated using p — A collisions 
with the y^-dependent cross-section mentoned above. In the present work, our 
motivation is to see the anomalous suppression due to Debye screening, we do not 
include nuclear suppression in our calculation. We also note here that the co- 
moving absorption will start when the plasma begins to hadronize. We have seen 
that the co-moving absorption is minimal with the parameters used in Ref. [27]. The 
extrapolation of these parameters to high energies is not straightforward. Therefore, 
in our calculation we do not include the hadronic absorptions. The expected pr 
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- - unsuppressed: (6-9) fm 
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,1,1, 



10 12 14 16 18 20 22 24 26 28 30 
P T (GeV) 



Figure 3: px distribution of J/ip in Pb — Pb collisions at y/s = 5500 GeV for 
Td = l.lT c ,x = 0.3 and t% = 0.5 fm/c. 



distributions of J /if} at LHC energies are shown in Fig. (3) for < b < 3 fm and 
6 < b < 9 fm centralities with and without Debye screening. For Td ~ 1.1 T c we 
see that the J/ips with p T < 15(10) GeV are suppressed for < b < 3(6 < b < 9) 
fm centrality. However, the other value of Td (see Table I) the J/ips will not be 
suppressed for the same value of the hard fraction x. We take Td = 1.5T C in Fig. (4). 
It is seen that suppression occurs below pt ~ 6 GeV. The value of the hard fraction 
x is not known at LHC energies at this stage. The survival probability decreases as 
x increases, since in that case, the initial temperature will be high (see eqs.(15) and 
(16)) as demonstrated in fig. (5). 

For ip' the Td values are small compared to that of J/ip and hence if>' will be 
suppressed substantially. This is shown in Fig. (6) with Td ~ 0.15T C at two different 
centralities. It shows that even a very high p? (~ 30 GeV) if)' is suppressed. We also 
plot the pt spectra of if)' for average values of Td at three centralities in Fig. (7) to 
show the centrality dependence of heavy resonance suppression. We observe that the 
suppression is small in going from most central to peripheral collisions as expected. 
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Figure 4: Same as Fig. (3) for T d = 1.5T C . 



10 p 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 i 
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1U 4 6 8 10 12 14 16 18 20 22 24 26 28 30 

P T (GeV) 

Figure 5: Same as Fig. (3) for Td = 1.1T C for different values of the hard fraction 
with impact parameter in the range < b < 3 fm. 
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P T (GeV) 

Figure 8: Same as Fig. (3) at y/H = 200 GeV, n = 0.6 fm/c, x = 0.1 and < b < 3 
fm. 

The prediction at RHIC energies is shown in Fig. (8). It is seen that J/ip is not 
suppressed even for most central collisions as the initial temperature is very close 
to the dissociation temperature. However, ip 1 is suppressed and the critical values 
of pt depend on the the dissociation temperature used. 

IV. Summary and discussions 

In this work we have calculated the pt distribution of charmonium by using colour 
evaporation model ( LO with gluon splitting) for RHIC and LHC energies. To val- 
idate our formalism we first reproduce the CDF data for px distribution of J/tp, ip', 
and Xc using CEM. It is assumed that at LHC energies QGP is formed and the heavy 
resonances are suppressed in a QGP due to Debye screening. We have calculated 
the survival probability using the prescription of Ref. [25] and convoluted it with 
the pt distribution of various resonances. First, the pt spectra of various charmo- 
nium states are obtained in nucleon-nucleon collision and then these are convoluted 
with Glauber model to get the same in nucleus-nucleus collisions. The centrality 
dependence of transverse momentum distribution is also estimated. We have seen 
that J ftp with pt > 15(10) GeV will not be suppressed in nuclear collisions for the 
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impact parameter in the range < b < 3 fm (6 < b < 9) fm. We do not include 
the co-moving suppression due to hadrons, as this is found to be small within the 
existing model parameters [27]. 
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